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1. Introduction 



The acceleration of the expansion of the Universe [1,2] has received no theoretical expla- 
nation yet. In these lecture notes, I will present an approach using scalar fields which can 
be motivated either from the dark energy [3] or the modified gravity points of view [4]. 
I have not attempted to cover the many aspects of the subject which can be found in 
very good review papers [5]. I have tried to emphasize the unity of the modified gravity 
models with a screening mechanism such as chameleons, dilatons and symmetrons, using 
as a template the example of large curvature f(R) gravity [6]. I have voluntarily excluded 
from the analysis the cases of the Galileon [7] and massive gravity [8] , with their associated 
Vainshtein screening mechanism [9]. This would have required a special treatment which 
can be found elsewhere [10]. 



2. Motivation 

The Hubble diagram of supernovae of the type IA first showed that the acceleration pa- 
rameter 

90 = T2", (2.1) 

a Q 

where ao is the scale factor of the Universe 1 , is negative implying that Hq > and the 
Universe accelerates. Assuming that the Cosmological principle stands and that General 
Relativity (GR) describes the large scale structure of the Universe, the combined data from 
Cosmic Microwave Background (CMB) observations, large scale structures and Baryon 
Acoustic Oscillations (BAO) lead to the energy contents of the Universe where Q m Q ~ 0.3 
and 17ao ~ 0.7. These two numbers are the energy fraction in matter and dark energy. 
Assuming that GR is valid, the Raychaudhury equation yields 

- = ^(P + 3p) (2.2) 

Ob O 

where Gat is Newton's constant, p the energy density of the Universe and p the pressure. 
Ordinary matter is characterised by an equation of state 

p = wp (2.3) 

where w = for non-relativistic matter and w = 1/3 for radiation. Obviously, ordinary 
matter with a positive equation of state cannot lead to the acceleration of the expansion of 
the Universe. Hence one must question one of the two hypotheses underlying this result. 
One may introduce a new type of energy such that the equation of state becomes negative 
w < —1/3: this new type of matter is called dark energy. On the other hand, if the 
acceleration of the Universe results from a misunderstanding of the laws of gravity on large 
scales, gravity must be modified and GR altered. 

In both cases, the acceleration of Universe occurs at very low energy of the order of 
the critical energy of the Universe p c ~ 10~ 48 (GeV) 4 implying it should be describable 



1 An index indicates the present Universe. 
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within the framework of low energy effective field theory. If this were not the case, the 
acceleration of the Universe would be at odds with all the rest of modern particle physics, 
which underlies the description of the Universe after Big Bang Nucleosynthesis (BBN). 
Within this context, an effective field theory should be parameterised by a few couplings 
and masses to be measured experimentally. Weinberg's theorem states that the unique low 
energy field theory of spin 2 particles respecting Lorentz invariance is GR [11]. As Lorentz 
invariance seems to be a fundamental fact of the Universe, the easiest way of modifying 
gravity at low energy is to introduce a new field beside the usual graviton. The simplest 
modification of gravity correspond to adding a mass to the graviton 2 , this automatically 
generates a field of scalar helicity. Similarly, taking cues from the acceleration of the early 
Universe, i.e. inflation, dark energy is most easily modeled out using scalar fields. Hence in 
both the dark energy and the modified gravity cases, scalar fields seem to be compulsory. 
We will describe how this happens in practice using our first example: f(R) gravity. 

3. f(R) Gravity 

If GR is not the theory of gravity on large scales, maybe a theory described by 

where f(R) = R + h(R) (and h(R) = for GR) is the right one on large scales? This 
theory is one amongst many like 

d XV ~ 9 16^ [6 - 2) 
which generically suffer from the Ostrograski instability whereby the Hamiltonian of the 
model is unbounded from below [12]. f(R) gravity does not suffer from this problem. This 
can be most easily seen by redefining the metric 

9,u = e W/mpI <7* (3.3) 

where the metric g^ v is the Einstein frame metric. The Einstein frame is always identified 
by choosing the metric allowing one to write the Lagrangian as 

By redefining the metric this way, as expected, a scalar field appears and the action becomes 

where the contractions are to be taken with the Einstein metric. In this frame, the theory is 
perfectly well-defined as long as the potential V(<j>) is bounded from below. Its expression 
is 

V(<t>) = (3.6) 

Z JR 



I 



2 The usual Pauli-Fierz theory of massive gravity is ridden with the Boulware-Deser instability in curved 
space, hence massive gravity is a lot more subtle and recent developments can be found in [8, 10]. 
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and the scalar field is obtained using the mapping 

fR = 1R = e_W/mP1 ( 3 - 7 ) 

where 8itGn = nipf and 

P-±. (3.8) 

A very useful example, which we will use as a template, consists of the large curvature 
models [6] where 

f(R) = R- W7rG N p A + (3.9) 
for R> Rq. Expanding to leading order we find that 



<A _f Ro Ro 



n+l 



mpi 2/3 

and 



(3.10) 



Notice that /?a plays the role of a cosmological constant and the potential is a decreasing 
function of (f> < mpi as long as i? > i?o- In the large curvature regime, this model provides 
an interesting example of modified gravity. We will investigate its properties in the course 
of these lectures. 

The acceleration of the Universe occurs when the potential V{4>) is a dark energy 
potential. Generically, dark energy models require that the mass of the scalar field <p is of 
the order of the Hubble rate 

m ~ H ~ 1(T 43 GeV. (3.12) 

This is a very low value implying that the range of the scale field is of the order of the 
present cosmological horizon. This would not be a problem at all if the scalar field (j) now 
were decoupled from matter. In most relevant cases, this is not the case in particular when 
the f(R) models are defined in the Jordan frame of matter 



S 



= j d^x^-g-^- + S m (^ 9lu/ ) (3.13) 



where S m is the standard model action for matter particles and the coupling of particles 

to gravity is mediated by the gravitons arising from the metric g^ v . In the Einstein frame, 
this matter action becomes 

S m (i>,9^) -> S m ^,e 2 ^ m ^g* u ). (3.14) 
This implies that fermions couple to the scalar fields as 

C = P-^-mMtpip (3.15) 
mpi 
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where m^o is the mass of the fermion in the Jordan frame. For a nearly massless scalar 
fields the propagator is — ^ £ and the tree level Feynman diagram with a scalar exchanged 
between two fermions gives a potential 

S^ N = -2? °"™*° (3.16) 
r 

which corrects Newton's potential by a factor 2/3 2 . Such a correction is tightly constrained 
by the Cassini measurement [13] 

/3 2 < 1(T 5 (3.17) 

which is violated by f(R) gravity. This does not mean that f(R) gravity is ruled out as 
the scalar force is screened in dense environments. 

4. Scalar-Tensor theories 

f(R) gravity in the Einstein frame is an example of a scalar-tensor theory defined by 

s = j ^yz^ ( _^|_ _ M! _ vm + Sm {^A\^ v ) (4.i) 

where A(<f>) is an arbitrary function. The dynamics of these theories are governed by the 
Einstein equation 

V - \^uR = 8wG N (T£ + Tp (4.2) 

where we have suppressed the E index for convenience. The scalar field has the energy- 
momentum tensor 

1% = d^dA - g^yyf- + V^)) (4.3) 
and the Klein-Gordon equation becomes 

= _ J^T m (4.4) 
where T m is the trace of the matter energy momentum tensor. We have defined the coupling 

h = m Pl ^ ' . (4.5) 
Matter is not conserved anymore as 

£> T^" = -^T m 9>. (4.6) 
mpi 

This can be better understood in the case of non-relativistic matter 

T% = peuW (4.7) 

where = ^ the velocity four- vector, r the proper time and u^u^ = — 1. The non- 
conservation of matter can be reexpressed as 

PE + Zh PE = ^p E <P (4.8) 
mpi 
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where the local Hubble rate is such that 3h = D^u^, and pE = u^D^pe- This is comple- 
mented with the Euler equation 

{d^ + 4>u^). (4.9) 



mpi 



We will come back to this relation later, let us first concentrate on the conservation of 
matter. It is particularly useful to define the conserved matter density p such that 

p E = A{4>)p (4.10) 
which satisfies the usual conservation equation for non-relativistic matter 

p + 3hp = 0. (4.11) 
The Klein-Gordon equation takes then the very convenient form 

D 2 (f) = ^ (4.12) 

where the effective potential is 

V cS (<f>) = V(<t>) + (A(<f>) - l)p. (4.13) 

With a decreasing V(</>) and an increasing A((f>), the effective potential acquires a matter 
dependent minimim (f) m \ n {p) where the mass is also matter dependent m(p). As the density 
increases, the mass increases in such a way that the effects of the scalar field is screened in 
a dense environment [14]. We will come back to this point later. 

Let us consider the large curvature f(R) models, in this case we find that the minimum 
is located at 

where pt = m^Ro = \p^ + Q, m opo and po = SH^m^y Notice that 4> m i n (p) is smaller 
is denser environments. When matter is the cosmological matter density which varies as 
p(a) = VL m Q^ and the redshift is 1 + z = ^, the mass at the minimum is simply given by 

m(p) = mo( -ttz — = ) l + )f (4.15) 

where 



/4^ao + $7 m0 , a s 
m = H 0X - — — — — (4.16) 
V (n + l)/ Ro 

Notice that the mass now is larger than the Hubble rate when is small enough. 
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5. Gravitational tests 

Scalar-tensor theories lead to a modification of gravity which can be easily analysed using 
the Euler equation in the quasi-static limit where \<j)\ <C \d(f)\. In this case we find that 

u, = -^d,4>. (5.1) 
mpi 

In the non-relativistic limit where u l = v l is the velocity field of the particles, and the 
metric in the Newton gauge is 

ds 2 = -{l + 2$ N )dt 2 + dx 2 {l -2$jv) (5.2) 

where <&n is Newton's potential, we find that the Euler equation reduces to Newton's law 



dv l 
~dt 

where the Newtonian potential is modified 



-0** (5.3) 



* = <Z> N + \nA((/)). (5.4) 

Hence gravity is modified by the presence of the scalar field. In the f(R) gravity case we 
have 

^ = <S> N + /3—. (5.5) 
mpi 

The same result can be obtained by writing the relativistic interval in the Jordan frame 

ds 2 j = A 2 (<p)ds 2 (5.6) 

as 

ds 2 = -{l + 2^)dt 2 + dx 2 (l-2$) (5.7) 

where 

$ = ^ N -\nA{4>). (5.8) 

Hence non-relativistic particles, which couple minimally to the Jordan frame metric feel 
the presence of two Newtonian potentials $ and ^ . This has important consequences as 
light bends according to $ + ^ = 2<3?tv which does not depend on the scalar field while 
matter particle evolve in the gradient of *S>. In cosmology, this implies that the velocity 
field of galaxies is sensitive to ^ and therefore to the scalar force while cosmic lensing only 
depends on <3?tv independently of the scalar field. 

Let us now analyse how gravity is modified for spherical and static objects. 

5.1 Point-like particle 

For a point-like particle, the matter density is simply 

p = m5^\r) (5.9) 
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where m is its mass. The Klein-Gordon equation for a theory like f(R) gravity where the 
coupling /3 is constant becomes 



implying that 



and 



d 2 (P = — mSWfr) (5.10) 
mpi 



Mr) = m (5.11) 

47rm P ir 



Hence we retrieve the fact that the Newtonian force is larger by a factor of (1 + 2/3 2 ). 
5.2 Small objects 

Let us consider a spherical object of density p c embedded in a sparse environment of density 
Poo. Far away from the object, the field converges to the minimum of the effective potential 
0oo in the density p^. Inside the object, the field is a small perturbation around a value 
4>o- At distances smaller than the range m^ 1 , the field outside reads 

<fr = 0oo + - (5.13) 
r 

and the field inside is 

Vn .sinhmor _,, , ril . 

<t> = 00 + -| - I) (5.14) 

rag mor 

where mo is the mass at (fio and V ' the derivative of the effective potential. The potential 
outside is defined by 

D=4(coshm ii-^^) (5.15) 



and 



— ^-(l — coshm i?) = (poo — (fro- (5. 16) 

m 



When the object is small and its radius R satisfies moR < 1, we find that 

00 ~ 0oo - 3Pm Pl $ N (R) (5.17) 

where $n(R) is Newton's potential at the surface of the body. As the object grows, (fro 
becomes smaller and smaller until it becomes almost uniformly equal to the minimum of 
the effective potential (fr c for the inside density p c . For the small radius case moR <C 1, the 
total Newtonian potential outside is simply 

* = (l + 2f3 2 )$ N (r) (5.18) 

which is the same result as in the the point-like particle case. 
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5.3 Thin shell 

When moR 3> 1, the field is almost uniformly equal to 4> c inside the object, apart from a 
thin shell at the surface of the object [14, 15]. The solution outside is 

( f> = <t> 0O -2ft eS m Pi <$> N (R)- (5.19) 

r 

where 

A /? 

fteS = 3/3— (5.20) 

and the thin shell extension satisfies 



AR 4>oo - (f> 



(5.21) 



R 6ft<f> N (R)m Pi ' 
The Newtonian potential outside the object is 

* = (l + 2ftft eff )<S> N (r) (5.22) 

implying that the modification of gravity is screened as long as 

|<Aoo - <Ao| < 2ft<S> N (R)m Ph (5.23) 

This is the screening criterion for models like f(R) gravity with a constant ft. In fact, 
the same criterion holds for all models with screening such as chameleons, dilatons and 
symmetrons where ft should be replaced by ftoo [16]. 

In practice, when one applies these results to the Cassini experiment, one must impose 

that 

/3oo/3cff < KT 5 (5.24) 

Another type of constraints can be obtained from cavity tests of the gravitational force. 

Inside a cavity of radius R filled with a vacuum density and surrounded by a bore 
of density p c , the field is essentially constant close to the centre of the cavity with a value 
0o determined by 

sinhmpfi _ o mg 
1 + ^R~-~^f (5 - 25) 

When this condition is satisfied, the test bodies inside the cavity must be screened and the 
effective coupling satisfy 

Pes < 1(T 5 (5.26) 

where 

2$ N (R)mpi 

and $at is the Newtonian potential of the test particles. 

In the case of large curvature f(R) gravity, the Cassini constraint reads 

f Ro < 10~ 5 $ JV (sun)(^)" +1 (5.28) 
Pt 
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where ~ W 6 pT and $7v(sun) ~ 10 6 implying that 

f Ro < 10 6n " 5 . (5.29) 

This is a very mild constraint. The cavity constraint is even milder as there is no solution 
for 0o implying that <p c is the value inside the cavity and therefore the test particles are 
effectively decoupled from the scalar field. 

6. Models 

Large curvature f(R) gravity is not the only type of models of interest. We will sketch 
three different ones here. 

6.1 Chameleons 

Chameleons [14, 15, 17] share the same coupling as f(R) models with 

A(4>) = eP^ (6.1) 

where (3 is a free parameter. Typically, the potential for chameleon models can be taken 
of the inverse power law form with 

A n+4 

= A 4 + — (6.2) 

where A ~ 10 -3 eV both acts like a cosmological constant and implies that the Cassini 
and cavity bounds are satisfied thanks to the thin shell effect. 

6.2 Dilatons 

Gravity tests are evaded in a very different way in the dilatonic models [18]. This is due 
to the shape of the coupling function 



Mo) - L + -^(0-^)2 (0.3) 



which implies that 

h = — {<t>-4>*) (6-4) 
mpi 

and the coupling to matter converges to zero when <fi is stabilised close to cf)± in a dense 

environment. The potential for dilaton models must be smooth, positive and slowly varying 
like 

V(<p) = V e-' f) / mpi (6.5) 

where Vo gives the order of magnitude of the energy density required to generate the 
acceleration of the Universe. The minimum of the effective potential 0min(p) — )■ 0* in a 
dense environment where (3(<f>*) = 0. 
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6.3 Symmetrons 

The symmetrons [19-21] have the same type of coupling as the dilatons where 

m = 1 + 2^ff (6 - 6) 



implying that 



mpi 

The potential is a Mexican hat with 



(6.7) 



V(4>) =V - ^<ft 2 + \<ft A (6.8) 
implying that there is a symmetry breaking transition for a density 

* = ^ (6.9) 

for which in an environment with p > p+, the minimum is at the origin with a vanishing 
coupling 

^min(p) = -> /3 (0 min ) = 0. (6.10) 

Hence gravity is only modified at low density, corresponding to the late time Universe or 
sparse environments. 



7. Cosmology 

The cosmology of screened models of modified gravity is universal at the background level 
and differs only for perturbations. 

7.1 Background Cosmology 

When the mass of the scalar field at the minimum of the effective potential is large enough, 
i.e. such that m(p) 3> H, the minimum is stable. This implies that if the field settles at 
the minimum early enough in the Universe, it will stay there for the rest of the evolution 
of the Universe. This is particularly important as the particle masses have a (ft dependence 

J7ty = A((ft)m^ (7.1) 

where is the mass in the Jordan frame. If A((f>) varied abruptly during Big Bang 
Nucleosynthesis (BBN), the formation of the elements would be altered. If the field follows 
the time evolution of the minimum before BBN, its stability is guaranteed especially when 
species like the electron decouple leading to a jump of the trace of the energy momentum 
tensor. If the field were not at the minimum, this would lead to kicks to the field values 
and therefore a large variation of the particle masses. This is not the case when the field 
follows the minimum from a redshift around z ~ 10 10 . 
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At the background level, the energy density of the scalar field is 



^ = y+Veff(0) (7.2) 



and the pressure 



p* = y - v(<fi). (7.3) 



The effective equation of state [16] 



l + i^ = l + ^«£ + (A(*)-l& (7.4) 



is extremely close to —1 as 

* 2 _„,n »2, « 4 



and 



f^m „9 / ^ \ fir 



i)^ M3nmfl J(_)^ (7.6) 
implying that 

which is tiny number, and therefore these models behave like A-CDM as long a is not 
too small. This is in particular the case in the late time Universe. 

7.2 Perturbations 

The only hope of distinguishing the screened modified gravity models from A-CDM resides 
in the very different properties of cosmological perturbations. Let us come back to the 
conservation equation which we will study in the Newtonian gauge 

ds 2 = a 2 (??)(-(l + 2<S> N )dr] 2 + dx 2 (l - 2Q N )). (7.8) 

Putting 

u° = a -\l - $jv), u* = a"V (7.9) 

and 

where 9 = div 1 and H = a' /a with ' = d/dr], we find that 

s' = -e (7.ii) 

with 5= 5 -f. 

Similarly the Euler equation becomes 

v'i + Uvi + rtdij = -d$ N - -^-di<f> (7.12) 

mpi 
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which we linearise to obtain 



9> + U6 = k 2 <5> N + ^k 2 5(t> (7.13) 
mpi 

in Fourier modes. Neglecting the energy density of the scalar field, the Poisson equation 
becomes 

**» = ->*3r- (714) 

In the quasi-static approximation, the Klein-Gordon equation becomes an algebraic relation 

5(f) _ f3<pp5 



Combining these equations we find that the density contrast must satisfy [17] 



(7.15) 



S " + - H s'-^^-(l + e(k,a))S = (7.16) 

where 

2/3 2 

e(k,a)= J (7.17) 

captures the effects of modified gravity at the linear level. It can be immediately inferred 
that on large scales k <C am(p), the evolution of the density contrast is the same as in the 
A-CDM paradigm, i.e. in the matter dominated era it follows 

6~ a (ri). (7.18) 

In the radiation era, the density contrast grows logarithmically still. The main difference 
occurs on small scales inside the Compton wavelength of the scalar field k 3> am(p). In 
this regime the growth is enhanced with a growing modes scaling like 

5 ~ a v l 2 (7.19) 

where 

-l + A /l + 24(l + /3j 



v = . (7.20) 

This results is of course only valid when (3^ is constant, on the other hand the anomalous 
growth is present for all models even when depends on the scale factor via the evolution 
of the field </> m in(p)- We will see in the next section that the mass scales are such that 
the main effects of modified gravity arise below 10 Mpc where non-linear effects must be 
taken into account. This is investigated via N-body simulations, showing a distorsion of the 
matter power spectrum at the percent level which may be within reach with forthcoming 
surveys such as EUCLID. 
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8. Modified Gravity Tomography 



The screened models that we have presented so far are all defined by a Lagrangian and two 
functions V{(fr) and A((fr). The path from this Lagrangian formulation to observations such 
as gravity tests or the growth of large scale structure is not direct. From a phenomenological 
point of view, it would be more efficient to define models using quantities which are closer to 
the physical observables. This is easily realised if the mass m(a) and the coupling constant 
/3(a) are given as a function of the scale factor of the Universe. Indeed these two functions 
completely characterise the time evolution of the linear cosmological perturbations. In 
particular, the mass m(a) indicates which scales are affected by the anomalous growth 
of structure due to the scalar field. In terms of matter density, the cosmological matter 
density p(a) varies from 10 g/cm 3 at redshifts around z = 10 10 before BBN to cosmological 
densities now for z = 0. Hence this parameterisation is in one to one correspondence with 
all the range of densities accessible to experiments in the solar system and in astrophysics 
3 . In fact the knowledge of m(a) and /3(a) is entirely enough to fully define the non linear 
functions V((f>) and fy. This is simply a reconstruction mapping which is a mathematical 
property, it has nothing to do with the fact that the field follows or not the minimum of 
the effective potential since before BBN. 

The field values can be parametrically defined by the following expression [22] 



where <p c is the value of the field at the minimum in a dense region of density p(a; n i). 
Taking a m i ~ 1CT 10 implies that this value of the field is the one inside dense bodies on 
earth where p(a- m i) ~ 10 g/cm 3 . The value of cfr(a) is the minimum one in a body of 
density p(a). Similarly the potential can be obtained using 



where Vq is a constant. Eliminating a from these two expressions, one obtains the potential 
V{4>) and the coupling (3^. This method is particularly important for numerical simulations 
as models defined by m(a) and /3(a) can be easily analysed in their non- linear regime using 
the reconstruction mapping. 

This method is also very efficient to impose the screening condition. Indeed it can be 
reformulated as 



where (fro is the value of the minimum far away from the object where the coupling is /3o 
and the object has a Newtonian potential at its surface equal to <3?g and a value at the 
minimum equal to (fro- This can be reexpressed a 




(8.1) 




(8.2) 



\(fr G - | < 2/3 $ G 



(8.3) 




(8.4) 



3 Larger densities like in neutron stars require to know the mass and the coupling for a much large z, 
this is not necessary to what follows. 
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where we have defined m(a) = mo/(o) and /3(a) = (3og{a). Let us apply this inequality to 
the Milky Way which must be screened to avoid a disruption of the dynamics of the galactic 
halo. In this case ac = 10~ 2 corresponding to a density inside the galaxy which is 10 6 larger 
than the cosmological density and ao = 1 assuming that the Milky Way is surrounded by 
the cosmological vacuum. For models like f(R) gravity in the large curvature regime where 
the integral da a /j"( a ) = and upon using that the galactic Newtonian potential is 

$ G ~ 1(T 6 , we get that [22,23] 

^>10 3 - (8.5) 

This condition is independent of f3o and means that any screened modified gravity model 
will have effects on Mpc scales only. For the case of large curvature f(R) gravity, this 
implies that [6] 

fn < 10- 6 (8.6) 

This is a loose bound as the screening of the Milky Way is not a strong quantitative 
constraint. A slightly stronger bound fn < 5 • 1CT 7 comes from distance indicators of 
screened vs unscreened astrophysical objects [24]. 

9. Quantum Corrections 
9.1 Effective field theory 

The modified gravity models with a screening property are at best effective field theories 
valid below the electron mass where all the massive fields of the standard model of particle 
physics have been integrated out 4 . This is an appropriate description of cosmology after 
BBN as the matter particles can be appropriately modeled using a fluid approximation. 
As a result the only quantum fluctuations which are present in this low energy model are 
the ones of the scalar field itself whose mass is extremely low below the cut off around the 
MeV scale. To evaluate these quantum corrections one must expand the field around the 
background value 4>o of the scalar field in the presence of matter. The effective potential 
reads then 

2 to 

where c p are dimensionless coupling constants 



(9.2) 



This expansion can also be written 

Wo + H) =V + ^fW 2 + A 3 ^ 3 + + £ (9.3) 

p>5 A P 

where each of the scales A p , p > 4 act as a cut-off scale for each of the higher order 
operators. For the non-renormalisable operators p > 4, these scales give an effective order 



4 Apart from the neutrinos. 
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of magnitude of the effective cut off of the theory as processes at energies larger than A p 
would violate the unitarity of the theory. The validity of perturbation theory also requires 
that |A| < 1. 

Let us consider a simple example where 

mpi (pT 

where e = — 1 when a > and vice versa. In the large curvature f(R) gravity case, we 
have the identification 

Vq = pa, Vi = 4p A , a = ^py, Ao = -^-JRo m P\ R 0i <PT = — • (9.5) 

For a constant coupling /3, we find that the minimum of the effective potential is located 
at 

*-*< e ^> 1/< - I> - « 9 ' 6 ' 

For these models we have 

A p ~ 4 = ^ (9 7) 

p 4!e(a-4)...(a-p+l) A 1 ' ' 

and 

|A| = (a-2)(a-3)^ (9.8) 

where 

ml = ea{a-l)^y-\ (9.9) 
As a result, the validity of perturbation theory is guaranteed when 

m < <P (9.10) 

and the effective cut-off is 

A p ~ <fo. (9.11) 

Therefore he validity of perturbation theory is not violated as long as the scalar field is light 
enough compared to the effective cut off. This guarantees the consistency of the model as 
the scalar field does not need to be integrated out to obtain the low energy effective theory 
below the cut off scale 4>o. We have also 

$-| (a -l)|*£t£) (9.12) 
and one can directly see that the theory is strongly coupled at high enough density. 
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9.2 Coleman- Weinberg corrections 

The one loop correction to the scalar potential at the one loop level is given by 

< 913 » 

where \i is a renormalisation scale which can be taken at the dark energy scale fi ~ 10~ 3 eV. 
This new potential may upset the properties of the screened modified gravity models [25] 
if its order of magnitude is larger than the dark energy scale now, if the minimum of the 
potential is shifted or even disappears and finally if the mass at the new minimum is much 
larger than m\. The first problem is easily solved as long as mo/Ho is not large, say 10 3 , 
as the Hubble rate now is so small. The order of magnitude of the correction to the first 
derivative of the potential at 0o is given by 

SV' ~ ^ ~ XVn. (9.14) 

As long as perturbative unitarity with |A| <C 1 is preserved, the shift of the slope of the 
potential is negligible. Similarly we have 

8ml ~ Xm o ( 9 - 15 ) 

and the mass shift is also tiny when perturbative unitarity is valid. Hence we find that as 
long as |A| <C 1, which is necessary to guarantee that the effective field theory makes any 
sense, the one loop corrections due to the scalar field are irrelevant. 

There is only one type of fermion which has not been integrated out, i.e. the neutrinos. 
The masses of neutrinos depend on the scalar field according to 

m-,p = A((f>)m^,o (9.16) 

implying that the one loop contribution to the scalar potential is 

64tt z mp\ 

This leads to a renormalisation of 

tfc-W,-^ (9.18) 

for each neutrino species. As < 10~ 3 eV in most probable scenarios of neutrino 

physics, we find that the effect of the neutrinos is to renormalise the cosmological constant 
by a value compatible with its observable value. 

10. Lorentz Violation 

We have focused so far on modified gravity where a scalar field couples to matter in a 
"conformal" way via the coupling function A((f>). This is not the only possibility by far and 
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another type of coupling includes a disformal term [26,27]. It involves the matter action 
Sm(^,9fMu) where matter couples to the metric 

9^u = A + ^ 4 ), (10.1) 

M is a suppression scale and the Einstein metric with which the Einstein-Hilbert term 
is written. 

Defining the energy momentum tensor as T^ v = -j= and expanding the action to 
linear order, we find that the scalar field couples derivatively to matter 

+ 1 AvCffi^J*' + J d l xsj^7jA 1 (4>)C m (:4',A 2 (4>)sf,„). (10.2) 

As soon as 9^ = d^ftd^^ does not vanish due to the presence of matter, Lorentz invariance 
is broken and a Lorentz violating coupling to the matter energy momentum tensor is present 
in the model. In a static configuration and in the presence of non-relativistic matter, the 
disformal coupling has no effect and therefore cannot be constrained by static tests of 
modified gravity. On the other hand, it can modify the propagation of fermions in dense 
environments. 

The action for massless fermions is 

Sf = ~ J d 4 x^ l -(^D^ - (D^)j^) (10.3) 
leading to the energy momentum tensor 

r£ = -\{^i {ix d v) ^ - (£(^ 7 ,)V>) (io.4) 



symmetrised over the indices. This induces the following interaction terms with the scalar 
field 

~2 d X ^~9E M 4 {^l( p D v) ip - (D^ip^ip). (10.5) 

In a typical situation where the metric is Minkowskian to a good approximation and the 
scalar field is static, the interaction term reduces to 



J d 4 x ^(i/ryidjif) - ditfrfjilf) (10.6) 



where d % = is a slowly varying function of space only. Hence a static configuration 
of the scalar field yields a Lorentz violating interaction in the Fermion Lagrangian. The 
resulting Dirac equation becomes 

i(-7°d + fdi + jcP-yidjW = 0. (10.7) 

The dispersion relation is obtained by squaring the modified Dirac operator to obtain 

Vl = (c 2 y a PlP a. (10.8) 
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This becomes the dispersion relation in an anisotropic medium with a square velocity tensor 

(c 2 ) ia = (5 ij + d?d j )(6 aj + d a dP). (10.9) 

The eigenmodes of the velocity tensor are d % and two vectors e\ , A = 1,2 orthogonal to d l . 
The eigenspeeds are Cd = (1 + \d\ 2 ) and twice c\ = 1. Hence, fermions go faster than light 
in the direction of the gradient d l (j), i.e. along the scalar lines of force, with 

Ac = c d - 1 = \d\ 2 . (10.10) 

This is of course something which has never been observed! This implies that the scale M 
must be large enough as the velocity of fermions could be larger than the speed of light 
in the thin shell of a screened object. This is very reminiscent of the wrongly announced 
result by the OPERA collaboration. . . 



11. Conclusion 



The models of screened modified gravity evade all gravitational tests. One may wonder how 
they could be efficiently probed. One possibility would be to find relevant situations where 
screening is not efficient. On astrophysical scales, this may happen in the Mpc range and 
may lead to a bump in the deviation of the power spectrum of density fluctuations from the 
A-CDM one [28,29]. The physics of stars, screened vs unscreened is also interesting [30]. 
Satellite tests of the equivalence principle may also be sensitive to an unscreened force [14] . 
In the laboratory, the Casimir effect experiments are extremely sensitive probes of new 
forces [31]. Finally, the interaction of slow neutral particles with matter, i.e. neutrons 
which are generically unscreened, may reveal surprises in low energy particle physics [32]. 
All in all, the search for modified gravity effects opens up a new era for physics at the low 
energy frontier. 
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